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Abstract 



We have searched for central production of a pair of photons with high 
transverse energies in pp collisions at yfs = 1.8 TeV using 70 pb" 1 of data 
collected with the D0 detector at the Fermilab Tevatron in 1994-1996. If 
they exist, virtual heavy pointlike Dirac monopoles could rescatter pairs of 
nearly real photons into this final state via a box diagram. We observe no 
excess of events above background, and set lower 95% C.L. limits of 610, 870, 
or 1580 GeV/c 2 on the mass of a spin 0, 1/2, or 1 Dirac monopole. 

PACS numbers: 12.60.-i, 14.80.-j, 13.85.Rm 



Submitted to Phys. Rev. Lett. 
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One of the open questions of particle physics is the existence of Dirac monopoles [|]||, 
hypothetical carriers of the magnetic charge proposed by P. Dirac to symmetrize the Maxwell 
equations and explain the quantization of electric charge. If such magnetic monopoles exist, 
then the elementary magnetic and electric charges (g and e) must be quantized according 
to the following formula: 



27m 

9 = , n = ±l,±2, 



(1) 



where n is an unknown integer. Here we assume that the elementary electric charge is that 
of an electron. If free quarks exist, Eq. (U) should be modified by replacing e with e/3, which 
effectively increases g by a factor of three. 

Dirac monopoles are expected to couple to photons with a coupling constant a g = 
g 2 /4:7T « 34 n 2 which is at least three orders of magnitude larger than the corresponding 
photon coupling to the electric charge (a e = e 2 /47r « 1/137). Therefore such monopoles 
could give rise to photon-photon rescattering via the box diagram shown in Fig. |I] ||f§]. 
The contribution of this diagram for pointlike monopoles to diphoton production at hadron 
colliders was recently calculated and shown to be significant even for monopole masses 
comparable to the collider beam energy. 

Since the virtuality (Q 2 ) of most incoming photons in the process of Fig. [I] is small |J, the 
interacting partons scatter at very small angles and therefore escape the detector through the 
beam pipe. Thus, a signature for monopoles at hadron colliders is the production of a pair 
of isolated photons with high transverse energies. This process gives a unique opportunity 
to find evidence for Dirac monopoles or to set limits on the monopole mass. 




FIG. 1. Feynman diagram for 77 production via a virtual monopole loop. 



The only previous study of this nature was made by the L3 experiment at the LEP 
e + e~ accelerator by searching for the Z — > 777 decay via a similar monopole loop J7|. It 
resulted in the lower 95% confidence level (C.L.) mass limit of 510 GeV for pointlike spin 
1/2 monopoles. Other accelerator experiments (see Ref. M) have looked for production 
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of monopoles by searching for the high ionization traces that would be produced by these 
particles, and would therefore be inherently restricted to monopole masses below the beam 
energy. A variety of experiments which look for monopoles in cosmic rays are sensitive to 
the relic monopole flux, rather than the monopole mass ||. Indirect limits on the Dirac 
monopole mass can be derived from measurements of the top quark mass and the axial and 
vector couplings of the Z to charged leptons ||. 

Despite numerous studies, QED with pointlike monopoles is still not a complete theory. 
For example, it is not clear whether such a theory can be constructed to be renormalizable 
to all orders ||. Also, arguments exist (see, e.g. 0) that Dirac monopoles must occupy a 
spatial volume of radius R ~ 0(g 2 /M), where M is the monopole mass, to accommodate 
the self energy implied by the large coupling. The non-observation of a new distance scale 
in QED or the SM for R < 0(1 TeV) requires the monopole mass to exceed ~ 100 TeV. 
Further theoretical work on this subject therefore is required to define the regions of validity 
for a theory of pointlike monopoles. In such a theory, it is possible that hard interactions 
of a monopole with photons would be weakened substantially by the effects of a monopole 
form factor. 

In this Letter we report on the results of a new search for Dirac monopoles with the 
D0 detector (described in detail elsewere |TIJ) operating at the Fermilab Tevatron proton- 
antipoton collider with beam energies of 900 GeV. The search is based on 69.5 ± 3.7 pb _1 of 
data recorded in 1994-1996 using a trigger which required the presence of an electromagnetic 
(EM) object with transverse energy Er above 40 GeV. This trigger did not require the 
presence of an inelastic collision, and therefore can be used to select low Q 2 events typical 
of the process in Fig. |I[ 

The following offline selection criteria are: (i) at least two photons with > 40 GeV 
and pseudorapidity |t/ 7 | < 1.1; (ii) missing transverse energy in the event $ T < 25 GeV; and 
(iii) no jets with E 3 T > 15 GeV and 1^1 < 2.5. The jet veto requirement is used to select 
the low Q 2 process in Fig. [I]. The trigger is > 98% efficient for this off-line selection. 

In order to determine the hard scattering vertex, we calculate the most probable direc- 
tion of each photon using the transverse and longitudinal segmentation of the EM calorime- 



ter |TTJ. These directions determine the position of the interaction vertex along the beam 
axis. The resolution on the vertex position for this method is 14 cm, taken from Z — > ee 
decays where the vertex can also be determined with high precision using the tracking in- 
formation. This EM-cluster-based vertex finding techique is preferred since for the event 
topology of Fig. [I] one does not expect charged particles, causing the tracking-based vertex 
finding to be biased significantly toward vertices from background interactions. We calculate 
kinematic parameters of the event based on the vertex obtained using the EM clusters. 



Each photon is required to have: (i) energy isolation |TTJ I < 0.1; (ii) more than 95% of 



the cluster energy deposited in the EM calorimeter; (iii) cluster shape consistent with that 
expected for a photon; and (iv) no tracks pointing toward the EM cluster from any of the 
event vertices. 

The overall efficiency for photon identification is (73.0 ± 1.2)% per photon, as detailed 
in Table |. This includes the (92 ±1)% probability of the photon not to convert in the 
material in front of the tracking chambers. The efficiency of criteria (i)-(iii) is determined 
using the Z — > ee events (with (ii) additionally checked using a geant |12[ simulation of the 
D0 detector for possible energy dependence); the efficiency of the no track requirement (iv) 
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TABLE I. Signal efficiency. 


uut 




rLtnciency (to J 




per photon 




r ^ n in 
1 < U.iU 




no n i n ^7 
»O.U ± U. ( 


EM fraction 




99.0 ± 1.0 


Shape consistency 




94.7 ±0.8 


No tracks 




91.1 ±0.4 


No photon conversions 




92.0 ± 1.0 




per event 




$ T < 25 GeV 




99.0 ±0.5 


Overall 




52.8 ± 1.4 



was determined using simulated photons obtained by rotating the electromagnetic clusters 
from Z — > ee decays by it/2 in azimuth [|TTJ. The overall efficiency for the diphoton selection 
is (52.8 ± 1.4)%. This includes the efficiency of the $ T veto (99.0 ± 0.5)% as well as the 
identification efficiency for a pair of photons. 

The above selection criteria define our base sample which contains 90 candidate events. 

The main backgrounds to photon scattering through a monopole loop are due to: (i) 
diagrams similar to Fig. [I] with other particles in the loop; (ii) QCD production of dijets 
(jj) and direct photons (jj) (with jets misidentified as photons due to fragmentation into 
a leading 7r° or rj decaying into a pair of spatially close photons, reconstructed as one EM 
cluster), or direct diphotons (77); and (iii) Drell-Yan dielectron production with electrons 
misidentified as photons due to tracking inefficiency. 



Background (i) is dominated by a virtual VT-loop and has been shown to be negligible [13 



The other two background contributions are estimated from the data. The QCD background 
is determined using the j'7 event sample collected with a single photon trigger, with the 
jet passing the same fiducial and kinematic cuts as the photon. We apply a a jet-faking 
probability P(j — > 7) which we measure to be (10.5 ± 1.5) x 10~ 4 by counting the number 
of photons in multijet events, and find the QCD background to be 25 ± 8 events. Direct 
photon and diphoton backgrounds are also included in this estimate. Their relative fractions 
are obtained from PYTHIA Monte Carlo (MC). The 30% error assigned to the QCD 
background estimate reflects the uncertainty in the direct photon fractions and in the jet- 
faking-photon probability. 

The Drell-Yan background is calculated from a sample of dielectron events passing the 
same fiducial and kinematic cuts as the signal sample. Multijet contamination of this sample 
is negligible since the probability for a jet to be misidentified as an electron is five times 
smaller than that for a photon. The probability for a dielectron pair to be misidentified as 
a diphoton pair is found to be (11 ± 1)% by comparing the number of events in the Z peak 
in the ee and 77 samples passing loose kinematic cuts. The Drell-Yan background in the 
base sample is 63 ± 7 events. The overall background in the base sample is 88 ± 11 (syst) 
events, in good agreement with the 90 observed candidates. 

To optimize the sensitivity of this search to the monopole loop contribution we apply a cut 
on the scalar sum of the transverse energies of all the photons in the event: St = Yli l±j 
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FIG. 2. Data and expected background as a function of S T nn cut. Points are data, the 
upper hatched region corresponds to the QCD background, and the lower shaded region shows the 
Drell-Yan background. The ~ 15% systematic error on the background is not shown. 



We vary the St cut threshold (S 1 ™ 111 ) in 10 GeV steps to achieve an expected background of 
0.4 events [16] . Such an optimization is based on the fact that for this expected background 



one has a 67% probability of observing no candidate events in the data in the absence of a 
signal. In such a case ||, the limits on the signal do not depend on the exact background 
value or its uncertainties. The agreement between the observed number of events and the 
predicted background as a function of S*™ 111 is illustrated in Fig. |2|. Note that since the plot 
shows the data and the backgrounds for St > S*™ 111 , the points are highly correlated. The 
S™ m = 250 GeV cut corresponds to a background of 0.41 ±0.11 events. The event in the 
base sample with highest St has St = 203 GeV, well below this cut. Taking into account 
the selection efficiency we set an upper limit for the production cross section of two or more 
photons with Y. e t> 250 GeV and |r/ 7 | < 1.1: 

a(pp ->> 27)| St >250 gcv,|^|<i.i < 83 fb (2) 

at the 95% C.L. This limit is obtained using a Bayesian approach with a flat prior and with 
the uncertainties in the efficiency and the integrated luminosity properly taken into account. 

Since the data are consistent with the background hypothesis, we can set limits on the 
production of pointlike Dirac monopoles. We calculate the acceptance for the monopole sig- 
nal using a fast MC program that generates diphoton events from a monopole loop according 
to the calculated differential cross section d? o j dE^drf 11 drf 12 || with a subsequent parametric 
simulation of the D0 detector. The MC model takes into account the interaction vertex dis- 
tribution; parton density distributions in the colliding protons and antiprotons, as described 



by the GRV [[17] parton distribution functions (p.d.f.); smearing of photon momenta; and 



detector acceptance. Figures |3](a) and (b) show the expected signal St distribution and the 
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correlation between the photon pseudorapidities, respectively. The cuts used in this analysis 
are indicated in the figures. The overall acceptance for the monopole signal is found to be 
(51 ± 1)%, where the error reflects variations due to different p.d.f. (estimated by taking 
the acceptance difference using GRV and CTEQ4L |1§), and uncertainty in the detector re- 



sponse parametrization. The acceptance does not depend on the monopole mass for masses 
above the typical photon energy (~ 300 GeV) |J. 

The total cross section of the heavy monopole production at the Tevatron is given by [Q : 

„( P p^ 77 + A-) = 57p(- g1 ^ J ) 8 fb, (3) 

where P is a spin dependent factor |[!,]13|]: P = 0.085, 1.39, and 159 for monopole spin of 
0, 1/2, and 1, respectively. The estimated error on this cross section due to choice of p.d.f. 
and to higher order QED effects is 30% |J. Additional uncertainties are associated with the 
77 — *■ 77 subprocess in Fig. [I] and with unitarity considerations. The coupling constant a g 
is replaced with an effective coupling || obtained by multiplying a g by a factor (_E' 7 /M) 2 , 
where E 1 is the photon energy, typically 300 GeV at the Tevatron. Both unitarity and 
perturbation theory assumptions are satisfied when this factor is <C 1 |3||| . 

Comparing the lower bound of the theoretical cross section corrected for acceptance with 
the cross section limit @ set by this measurement, we obtain the following lower limits on 
the pointlike Dirac monopole mass (see Fig. ^): 

( 610 GeV for S = 
M/n > I 870 GeV for 5 = 1/2 . 
{ 1580 GeV for S = 1 

These are currently the most stringent limits on heavy pointlike monopole mass. (We do 
assume, if more than one type of Dirac monopole exists, that there is no cancellation among 
the loop diagrams involving each monopole type.) 

We note that the effective coupling exceeds 1 and unitarity is violated close to the exper- 
imental bound. For values E 1 /M > 1, the cross section will grow more slowly, approaching 
the usual 1/M 2 behavior of a QED process || which satisfies unitarity. Also, for lower 
monopole masses the effective parameter of the perturbation theory used in the calcula- 
tions H becomes too large, and therefore one would expect a non-negligible contribution of 
the higher order diagrams with four, six, etc. photons in the final state. The latter effect 
is, however, largely compensated by our analysis cut on the sum of the photon transverse 
energies; if part of the signal cross section is due to the higher order diagrams, the above 
limits are unaffected. 

When more complete theoretical calculations are available, limits on the monopole mass 
could be updated by comparing the modified cross section expression with the experimental 
limit (0). With current theory || the above limits are strictly valid only for monopole masses 
above several hundred GeV. 

As a cross-check of the results of this search we have selected elastic or nearly elastic 
collisions by requiring no hits in the forward scintillating hodoscopes used for luminosity 



monitoring and triggering on the inelastic collisions [10]. This requirement drastically re- 



duces the backgrounds. The remaining background in the base sample for elastic events 
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0.5 1 -2-1012 

S T (TeV) r| Y/ 

FIG. 3. a) Normalized St spectrum and b) photon pseudorapidities for the diphoton produc- 
tion via a heavy monopole loop. The arrow in a) and square in b) show the chosen cuts in the 
corresponding parameters. 

is 1.8 ± 0.4 events, dominated by diffractive Drell-Yan events and residual inelastic back- 
ground due to inefficiency of the forward hodoscopes. We observe one candidate event in 
the base sample, consistent with this expected background rate. For S"™ 111 ~ 100 GeV the 
background is 0.4 events, and no candidates are observed. We use this sample only as a 
cross check because the efficiency of these selection requirements is significantly lower than 
that of the main analysis method, primarily because of multiple interactions. 

In conclusion, we have performed a search for heavy pointlike Dirac monopoles by search- 
ing for pairs of photons with high transverse energies. Our data agree with the expected 
background from QCD and Drell-Yan production. No candidates pass the final cuts. Using 
theoretical calculations B we set 95% C.L. lower limits on the Dirac monopole mass for min- 
imum magnetic charge (n = 1) in the range 610 to 1580 GeV, depending on the monopole 
spin. These are the most stringent mass limits on heavy pointlike Dirac monopoles to date. 
Our cross section limit (0) is 83 fb, and may be applicable to the other production processes, 
such as that of dyons M or other exotic objects strongly coupled to photons. 

We are grateful to I. Ginzburg and A. Schiller for many discussions and detailed cross 
section information and to U. Baur, B. Dobrescu, and A.S. Goldhaber for helpful discus- 
sions. We thank the staffs at Fermilab and collaborating institutions for their contributions 
to this work, and acknowledge support from the Department of Energy and National Sci- 
ence Foundation (U.S.A.), Commissariat a L'Energie Atomique (France), State Committee 
for Science and Technology and Ministry for Atomic Energy (Russia), CAPES and CNPq 
(Brazil), Departments of Atomic Energy and Science and Education (India), Colciencias 
(Colombia), CONACyT (Mexico), Ministry of Education and KOSEF (Korea), and CON- 
ICET and UBACyT (Argentina). 
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FIG. 4. The curved bands show the low and upper bounds on theoretical cross sections || for 
monopole spin, S = 0, 1/2, and 1. The horizontal line shows the 95% CL experimental upper limit 
(|2|) on the cross section. The arrows indicate the lower 95% CL limits on the monopole mass at 
each spin value. 
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